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A method for predicting thermodynamic effects of cavitation (changes 
in cavity pressure relative to stream vapor pressure) is presented. The 
prediction method accounts for changes in liquid, liquid temperature, 
flow velocity, and body scale. Both theoretical and experimental studies 
used in formulating the method are discussed. The prediction method 
provided good agreement between predicted and experimental results 
for geometrically scaled venturis handling four different liquids of widely 
diverse physical properties. Use of the method requires geometric simi- 
larity of the body and cavitated region and a known reference cavity- 
pressure depression at one operating condition. 


Because cavitation is a vaporization process that involves heat and 
mass transfer, the physical properties of the liquid and its vapor and the 
flow conditions can affect the cavitation process and thus the cavitation 
performance of hydraulic equipment as well. These combined effects of 
fluid properties, flow conditions, and heat transfer — termed thermo- 
dynamic effects of cavitation — can improve cavitation performance. For 
example, cavitation studies with pumps and pump inducers (refs. 1 
through 6) have shown that, for certain liquids and/or liquid tempera- 
ture, the net positive suction head (NPSH) requirements can be signifi- 
cantly less than those for room-temperature water. This improvement in 
inlet pressure requirements is attributed to the varying degrees of evapora- 
tive cooling associated with the cavitation process. Because of cooling, 
the vapor pressure of the liquid adjacent to the cavity and the cavity 
pressure itself are decreased relative to the vapor pressure of the bulk 
liquid. This decrease in cavity pressure retards the rate of further vapor 
formation, thereby allowing lower inlet pressures than are otherwise 
possible. For a given cavity size at fixed flow conditions, a reduction in 
cavity pressure allows a corresponding reduction in the inlet pressure 
requirements. The accurate prediction of the thermodynamic effects of 
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cavitation is therefore essential to an optimum system designed to 
operate with cavitation. 

The objective of this paper is to summarize the results obtained from 
theoretical and supporting experimental venturi cavitation studies (refs. 
7 through 11) which led to the formulation of a method for predicting 
thermodynamic effects of cavitation. The method was motivated by, 
but is an extension of, the 5-factor concept of Stahl and Stepanoff (ref. 
12) . The prediction method accounts for changes in liquid, temperature, 
flow velocit}^ and scale. The cavitation data were obtained for Freon-114, 
liquid hydrogen, liquid nitrogen, and water using two accurately scaled 
venturis (1.0 and 0.7 scale) as test devices. The liquid hydrogen data 
were obtained from a NASA-sponsored study conducted at the National 
Bureau of Standards (ref. 11). The data for the other liquids (refs. 7 
through 10) were obtained at the NABA Lewis Research Center. 

ANALYTICAL REVIEW 

Thermodynamic Effects of Cavitation 

Heat is required to generate the vapor within a cavitated region, and 
this heat of vaporization must be drawn from the surrounding liquid. If 
the time available for vaporization is short, as is usually the case for 
flowing systems, the liquid which supplies the heat for vaporization 
consists of a relatively thin film of liquid adjacent to the cavity. Cooling 
of this film causes the vapor pressure of both the liquid film and the 
vaporous cavity to drop from that of the bulk liquid to that corresponding 
to the reduced local temperature. For a cavity of given size with fixed 
flow conditions, a reduction in cavity pressure allows a corresponding 
reduction in the inlet pressure requirements. Thus the inlet pressure 
requirement for a particular flow device handling liquids which exhibit 
significant drops in cavity pressure is less than that for liquids which do 
not exhibit this beneficial effect. In general, cavity-pressure and tempera- 
ture depressions (thermodynamic effects of cavitation) are a function of 
fluid properties as well as the body geometry, the flow velocity, and the 
complex heat and mass transfer mechanisms involved. 

Heat Balance 

The magnitude of cavity-pressuni depression can be estimated by 
sitting up (as an initial step) a heat balance between the heat required 
for vaporization and the heat drawn from the liquid surrounding the 
cavity; that is, 

PrVrL = pt»tCl{AT) (1) 

where vi represents the volume of the thin layer of liquid that is cooled 
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during vaporization. Although applied to cavitation in a flowing system 
herein, the energy balance (eq. (1)) is equivalent to a static model of a 
simple insulated cylinder-piston model (ref. 13) . A given mass of liquid 
fills the cylinder. With an increase in total volume, achieved by lifting the 
piston, some of the liquid is vaporized. The heat required to form the 
vapor that fills the resulting increased volume must be drawn from the 
remaining liquid, thereby lowering the temperature by some amount AT. 
The primary assumptions associated with the heat balance (eq. (1) ) are 
discussed in reference 8. 

By expressing the temperature drop AT (eq. (1)) in terms of a vapor- 
pressure drop AK, through use of the vapor-pressure temperature rela- 
tions, equation (1) may be expressed in the form 


Pv Tj / dhv\ 

pi VI Cl \dT ) 


( 2 ) 


where dhy/dT is the slope of the vapor-pressure/temperature curve at the 
temperature of interest. All the terms in equation (2) are fluid properties 
and are known except the vapor- to liquid-volume ratio vy/vi. In the 
experimental case the absolute value of volume ratio cannot be directly 
determined because vi is only that volume of liquid actually involved in 
the vaporization process. Thus vi is only a fraction of the entire liquid 
stream. However, as shown in equation (2), the parameter Ahy can be 
expressed as a function of vy/vi for any liquid whose fluid properties are 
known. To account for changes in fluid properties as the temperature 
decreases due to vaporization, equation (2) may be solved (in terms of 
Ahv versus vy/vi) in incremental form. That is, suitably small increments 
of Ahv are assumed, and the results from successive steps are summed to 
yield Vy/vi. Results obtained in this manner for water, liquid nitrogen, 
Freon-1 14, and liquid hydrogen at various initial fluid temperatures are 
presented in figure 1. The nonlinearity of the curves is due to the change 
in fluid properties as the equilibrium temperature drops due to vaporiza- 
tion. 

For a given value of Vv/vi, the depression in cavity pressure can vary by 
several orders of magnitude depending on liquid and liquid temperature. 

By using the Clasius-Clapeyron equation to approximate the slope of 
the vapor-pressure/temperature curve, the following useful approxi- 
mation of equation (2) — and thus the curves of figure 1 — is obtained: 



This expression provides a useful and close approximation of the curves 
presented in figure 1. 
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a. — Water, liquid nitrogen, and Freon-114- 


b . — Liquid hydrogen. 


Figure 1. — Vapor pressure depressions for 
various liquids and liquid temperatures. 


temperature, yqijid nit,-Qgef, 

o' Water 



As previously discussed, the vapor-to-liquid-volume ratio vv/vi in any 
real flow situation is not known and cannot be measured directly. Thus 
values of the Py/vi ratio from figure 1 are used only in a relative sense. 
However, useful predictions for geometrically similar cavities may be 
made as follows. By experiment, an elTective or reference value of Vv/vi 
is established through the determination of cavity-pressure depressions 
Ah„ for one flow devir j, liquid, temperature, and velocity. Then values of 
Vv/pt for other geometrically similar flow devices, liquids, liquid tempera- 
tures, and flow velocities may bi; estimated relative to this reference value, 
as described subsequently. With these predicted values of p„/pi and figure 
1, determination of Ah„ relative to reference data is possible. 

Estimation of Vapor-to- Liquid Volume Ratio 

In the overall systematic study of thermodynamic effects of cavitation, 
several theoretical analyses for estimating vapor-to-liquid-volume ratio 
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were made. However, in the interest of brevity, only the pertinent .rela- 
tions and a short discussion of these analyses is presented. Details of the 
various derivations are presented in the references cited. 

The theoretical analysis of reference 8 shows that for one body scale, 
the vapor-to-liquid-volume ratio r »/ vi can be predicted for changes in the 
liquid used, the temperature, the velocity, and the cavity length, pro- 
vided that a reference value of v,/ vi is known. Subsequently this study was 
extended to include the effects of geometric scaling (ref. 9) . The combined 
results from references 8 and 9 led to the following general equation (see 
ref. 10) for predicting the vapor-to-liquid-volume ratio, relative to a 
reference value, for changes in liquid, liquid temperature, flow velocity, 
and model scale: 


M /AV"” [(AAl'T (4) 

V a / \Vo.rJ \DrJ [ Ax/D J 

where a is thermal diffusivitjq Fo is free-stream velocity, D is diameter, 
and Ax is the axial length of the cavitated region. The exponents m, n, 
and p depend, at least in part, on the heat transfer process involved and 
must be determined by experiment. The derivation and use of equation 
(4) requires, first, that geometric similarity of cavitating flow is main- 
tained for various liquids and flow conditions and, second, that thermo- 
dynamic equilibrium conditions exist within the cavitated region. At 
present, these two requirements can only be assumed for untested liquids. 

EXPERIMENTAL VENTURI STUDIES 

In order to determine the exponents for the theoretically derived 
equation (4), it was necessary to conduct cavitation studies in such a 
way that the temperatures and pressures within cavitated regions of 
various liquids could be measured directly. Accurately scaled venturis 
(1.0 and 0.7 scale) were used as test flow devices. In order to provide a 
wide range of physical properties, Freon-114, liquid hydrogen, liquid 
nitrogen, and water were used as test fluids. The fluids were studied over a 
range of approach flow velocities (100 to >200 ft/sec for liquid hydrogen 
and ~19-50 ft/ sec for the other liquids) . 

Venturis 

The venturis, made of transparent acrylic plastic, incorporated a 
circular arc contour to provide convergence from the approach section 
to a constant-diameter throat section (fig. 2). The throat section is 
followed by a conical diffuser. The venturi with the 1.743-inch approach- 
section diameter is referred to herein as the 1.0-scale venturi, and that 
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Figure 2. — Venturi test sections (all dimen- 
sions in inches). 
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with the 1.232-inch approach section diameter is referred to as the 0.7- 
scale venturi. Detailed descriptions of the venturis are given in references 
7 through 9. The venturis were mounted either in a small, closed-return 
hydrodynamic tunnel of 10-gallon capacity (NASA facility, ref. 7) or in 
an NBS liquid hydrogim blowdown test facility (ref. 11) with the 0.7-scale 
venturi located between 265-gallon supply and receiver dewars. 

The noncavitating wall-pressure distribution for the venturis is pre- 
sented in figure 3. As stream pressure was lowered, cavitation occurred 
first at the point of minimum pressure located on the circular arc portion 
of the venturi. The cavity’s leading edge was visually observed to remain 
fixed at the minimum pressure location for all cavitating conditions 
studied. 


Figure 3. — Wall pressure distribution 
for venturis. 
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Appearance of Developed Cavitation 

Photographs showing typical developed cavitation for Freon-114 and 
liquid hydrogen are presented in figure 4. Although the vapor formation 
in hydrogen was more homogeneous, the appearance of cavitation for 
these fluids was similar; that is, the cavities exhibit a uniform peripheral 
development with the leading edge at a fixed location. Because of this 
similarity in cavitating flow (a prerequisite in the analysis), the data 
obtained for both hydrogen and Freon-114 are considered equally valid 
for determining the exponents in the theoretical prediction equation. 
Temperatures and pressures within the relatively thin annulus of vapor 
formed (estimated thickness less than 0.03 inch in the throat) were 
measured simultaneously by pressure taps and calibrated thermocouples 
located on the adjacent wall. 


RESULTS AND DISCUSSION 

Specific trends in temperature and pressure depressions within cavitated 
regions of Freon-114 and liquid hydrogen are discussed first. Exponents 
in the theoretical prediction equation are then determined, and, finally, 
the significance of a developed cavitation similarity parameter Kc, min 
is discussed. 

Typical Cavity Pressure and Temperature Distribution 

A plot showing typical variations in pressure and temperature depres- 
sions within a constant-length cavity (Freon-114) is show in figure 5. 
Cavity-pressure depression is presented as a function of xjD measured 
from the minimum pressure location determined with noncavitating flow. 
The data are for a cavity length of about 2.75 inches (comparable to 
that shown in fig. 4d) and a constant approach velocity of 31.8 feet per 
second. All pressures within the cavitated region are less than vapor 
pressure based on the free-stream temperature (represented by = 0) . 
The maximum measured cavity-pressure depression always occurred at 
or near the cavity’s leading edge. The point at which the cavity pressure 
equals the vapor pressure of the bulk liquid (A/i„ = 0) coincides with the 
collapse region of cavitation. The open symbols represent cavity pressures 
measured directly, whereas the solid symbols represent the fluid vapor 
pressure corresponding to the locally measured temperature. The figures 
adjacent to the solid symbols represent measured reductions in local 
cavity temperature below the free-stream temperature. The good agree- 
ment between the measured pressures and the vapor pressures based on 
measured temperature indicate that, for engineering purposes, thermo- 
dynamic equibrilium conditions can be assumed to exist within the cavity. 
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a . — Liquid hydrogen in 0.7-scale ven- 
turi; nominal cavity length, 1 .0 free- 
stream diameter. 



b . — Liquid hydrogen in 0.7-scale ven- 
turi; nominal cavity length, 1.6 free- 
stream diameters. 


c. — Freon-114 in 0.7-scale venturi; 
nominal cavity length, I.4 free- 
stream diameters. 


d- — Freon-114 in 1.0-scale venturi; 
nominal cavity length, 1.5 free- 
stream diameters. 


Figure 4, — Typical cavitation in venturis. 
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Open symljols denote measured pressure 
depression 

Solid symbols denote vapor pressure 
depression at measured temper- 
ature depression 





Axial distance from minimum pressure location, x, in. 


Figure 5. — Typical pressure and 
temperature depressions within cav- 
itated region. Freon-114; 1.0-scale 
venturi; velocity, 31.8 feet per sec- 
ond; temperature, 519.1° R. 


Also, the studies showed that for fixed flow conditions and geometrically 
similar cavities, a change in the maximum cavity depression resulted in 
an equal change in the inlet pressure requirement (NPSH). This occurred 
for all liquids and temperatures studied. 


Effects of Temperature, Flow Velocity, Cavity Length, and Scale 

An increase in bulk liquid temperature caused an increase in cavity- 
pressure depression along the entire cavity length as illustrated in figure 6. 
For the conditions showm, the maximum cavity depression is about 8.5 
feet of Freon-114 (~5.5 psi) at 538.8° R, which is about 1.7 times that 
at 519° R and 4.6 times that at 467.5° R. This trend with temperature 
was observed for all test liquids at each flow velocity studied. However, 
for water, over the 500° R to 580° R temperature range studied, the 
cavity-pressure drop was so small as to be negligible. 

An increase in free-stream velocity caused an increase in cavity- 
pressure depressions. Some typical results are presented in figure 7. As 
velocity was increased from 22.9 to 44.5 feet per second, the maximum 
cavity depression increased from about 6.6 feet of Freon-114 to about 
11.0 feet. 

As cavity length was increased, the cavity-pressure depression in- 
creased over the full axial length of the cavity, as illustrated by the liquid 
hydrogen data (0.7-scale venturi) presented in figure 8. The trend shown 
is typical for the other liquids tested. For the conditions studied, the 
measured cavity-pressure depressions for hydrogen w'ere about 50 to 100 
times greater (in feet of liquid) than those for Freon-114. Maximum 
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Figure 6. — Effect of free-stream 
liquid temperature on cavity pres- 
sure and temperature depres- 
sions. Freon-1 H, 1.0-scale ven- 
turi; nominal cavity length, 2.75- 
inches. 
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Figure 7. — Effect of free-stream velocity on 
pressure and temperature depressions 
within cavitated regions of Freon-114. Ven- 
turi scale, 0.7; liquid temperature, 540° R; 
nominal cavity length, 1.6 inches. 
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depressions reached as much as 580 feet of hydrogen at the higher free- 
stream velocities and temperatures studied, compared to about 11 feet 
for Freon-114. 

For geometrically similar cavities (Aa:/D = const), cavity-pressure 
depressions over the entire length of the cavity are little affected by 
model scale, at least for the modest range tested. This is illustrated in 
figure 9, where P>eon-114 data obtained in both the 0.7- and 1.0-scale 
venturis are compared. Cavity-pressure depression is plotted as a function 
of the ratio of axial distance to free-stream diameter. 
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Figure 8. — Effect of cavity length 
on cavity pressure depression. 
Liquid hydrogen; 0.7-scale 
venturi. 
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Figure 9. — Effect of venturi scale on pres- 
sure and temperature depressions within 
cavitated regions of Freonrll4. Cavity 
length, 1 .4 free-stream diameters. 


DETERMINATION OF EXPONENTS IN PREDICTION 
EQUATION 

In all instances the maximum measured cavity-pressure depressions 
were used to determine the exponents for the equation for prediction of 
vapor-to-liquid-volume ratio v„/vi (eq. (4)). The Freon-114 data (for 
both the 1.0- and 0.7-scale venturis) and the hydrogen data (0.7-scale 
venturi) were used in conjunction with the method of least squares to 
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solve for all the exponents simultaneously. More than 100 data points 
were used. The resulting expression with the experimentally determined 
exponents is 



■ At/D 
_ (At/D) ref_ 


( 5 ) 


The exponent for the diameter-ratio term presently comes via theory 
from the exponent n of the velocity-ratio term (see eq. (4)). Thus addi- 
tional cavitation studies which use other scale factors are needed to 
further substantiate the exponent for D/D^et- 

From a measured value of A/i„ in a model of given scale, the volume 
ratio {v^/ vi) ref can be determined from figure 1 (or eq. (2) ) . Equation (5) 
can then be used to determine {vv/vi)pred for the liquid, temperature, 
velocity, and scale of interest. With this value of (r„/r;)pred, figure 1 can 
be used again to predict the corresponding value of A/i„. For the present, 
equation (5), with the experimentally determined exponents, is offered 
as the best approach toward generalizing results for different scaled 
venturis, different liquids and liquid temperatures, and various velocities. 


COMPARISON OF PREDICTED AND EXPERIMENTAL 
RESULTS 

The agreement to be expected between experimental and predicted 
results is presented in figure 10. The complete test range of Ali„ values is 
presented in figure 10a, which emphasizes liquid-hydrogen data. The low 
range of values is presented in figure 10b. All points shown in figure 10 
are based on a single reference point represented by the solid symbol. The 
reference value was obtained with 538° R Freon-114 in the 1.0-scale 
venturi at a free-stream velocity of 19.6 feet per second. The cavity 
length, At/D, was 0.72. The measured Ahy was 5.0 feet of Freon-114, 
and thus the value of (vv/vOrei (from fig. 1) was 2.6. Predicted values of 
Vv/vi were then calculated from equation (5). All predicted A/i„-values 
appear to be quite good, considering that predicted values of up to 600 
feet of liquid (at Fo>200 ft/sec), were based on a single reference value 
of only 5.0 foot at 19.6 feet per second. 

It is concluded that the simple static model analysis with experimentally 
determined exponents provides reasonably good predictions over a wide 
range of fluid properties and flow conditions, provided that (1) the 
maximum cavity-pressure depression is known for at least one operating 
condition, (2) similarity of cavitating flow is maintained, and (3) local 
equilibrium of temperature and pressure exists at the vapor/liquid 
interface. 
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Figure 10 . — Comparisons of measured and predicted 
vapor-pressure depressions. 


CAVITATION SIMILARITY PARAMETERS 


The conventional cavitation parameter is usually expressed as 


X.,= 


ho 

ToV2g 


( 6 ) 


The parameter for developed cavitation is derived from the assump- 
tion that Bernoulli’s equation for steady ideal flow applies between a 
free-stream location and the cavity surface and that the cavity surface 
is at a constant pressure equal to free-stream vapor pressure. However, 
because of the previously discussed thermodynamic effects of cavitation, 
the cavity pressure can be significantly less than free-stream vapor pres- 
sure and can also vary with axial distance. Thus, a more general expres- 
sion for the cavitation parameter would have in equation (6) replaced 
with a more appropriate reference pressure in the cavity. As in references 
8 through 10, the maximum cavity-pressure depression is selected as the 
reference pressure to define the following developed cavitation parameter: 
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ho /ic.min ^0 hv ^ (A/lv)max l (AAt»)max 

A„„i„= F„V2g ^ ’’■^>„V2g 

The venturi cavitation studies show that Kc, min is approximately con- 
stant over the range of liquids, temperatures, velocities, and venturi 
scales tested, provided that the geometric similarity of the cavitated 
region is maintained. The value of Kv however, varied widely. 

The degree to which a constant value can be used to represent Kc, min 
is shown graphically in figure 11, where the numerator of Kc,m\n is plotted 
as a function of its denominator. The line shown is for a K^.min of 2.47. 
The good agreement between the data and the line indicates that a single 
Kc.min value of 2.47 can be utilized to represent all the hydrogen and 
Freon-114 data. This Kc, min value of 2.47 also applies quite well for liquid 
nitrogen and water. As previously stated, the cavity-pressure depressions 
for water were practically negligible; thus, the values of Kc.min and the 
conventional cavitation parameter Kv were essentially the same for the 
range of water temperatures studied. Although the value of Kc.min re- 
mained constant for different cavity lengths, it is possible that for other 
venturi shapes Kc.min will vary with cavity length. 

Prediction of Free-Stream Pressure Requirements 

For constant values of Kc.min and Fo in equation (7), a change in 
(A^t.)max results in a corresponding change in ha— hr. This pressure 


Figure 11 . — Effects of venturi scale and 
liquid on minimum pressure of cavitated 
region. Various cavity lengths, Ax/I). 
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difference ho — hv is a measure of free-stream pressure requirements to 
obtain geometrically similar cavities and includes the thermodynamic 
effects of cavitation. Thus, for the wide range of conditions for which 
Kc,min is constant, it is a most useful cavitation similarity parameter. 
When at least one value of (AA„)max (or Ac, mm) is known and used as the 
reference value for predicting other values of (AA„)max, as previously 
discussed, the free-stream pressure requirements ho for untested conditions 
can be predicted by equation (7) . This equation can be rearranged so that 

Vo^ 

ho=-:^(K c.min) ^c,min ) (8) 


where hv he, min — 

In the present study the reference value of AA„, needed for prediction 
purposes, was measured directly. However, because Kc, min remains con- 
stant with geometric similarity of body and cavitated region, it is possible 
to estimate (AA,,)m-x=ref without measuring cavity pressure directly. In a 
given flow device, two different tests are made in which the free-stream 
pressure requirements {ho—K) to obtain geometrically similar cavitated 
regions are measured. These tests need not be for the same liquid, tem- 
perature, velocity, or model scale; however, at least one of these tests 
must yield measurable thermodynamic effects of cavitation (AA„>0). 
From these measured free-stream pressure requirements, it is possible to 
determine a (AA„)ref value using an iterative process as described in 
reference 14. 


CONCLUDING REMARKS 

A method for predicting thermodynamic effects of developed cavitation 
in venturis, or other flow devices, has been developed. The prediction 
method, formulated from simple theoretical analyses and supporting 
experimental cavitation studies, provided good agreement between pre- 
dicted and experimental results for geometrically scaled venturis handling 
liquids of widely diverse physical properties. The prediction method ac- 
counts for changes in liquid, liquid temperature, flow velocity, cavity 
length, and model scale, provided that similarity of the cavitating flow is 
maintained. 

The method requires a known reference value of the maximum cavity- 
pressure depression at one operating condition, obtained either by direct 
measurement or from two appropriate tests wherein stream static pressure 
is measured under conditions of similarity of the cavitating flow. These 
latter test data need not necessarily be for the same liquid, temperature, 
velocity, or model scale. 
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Although the present study is limited to venturis, the basic prediction 
method has been applied successfully to the prediction of cavitation per- 
formance for devices such as pumps and inducers (see ref. 14). 

LIST OF SYMBOLS 

Cp noncavitating wall-pressure coefficient, (hx—ho)/ {Vii‘/2g) 

Cl specific heat of liquid, Btu/(lb mass) (°R) 

D free-stream diameter (approach section, fig. 2) , in. 

g acceleration due to gravity, 32.2 ft/sec- 

hc, min minimum measured pressure in cavitated region, ft of liquid abs 
hv vapor pressure corresponding to free-stream or bulk liquid 
temperature, ft of liquid abs 

Ahv decrease in vapor pressure corresponding to decrease in liquid 
temperature, ft of liquid 
hx static pressure at x, ft of liquid abs 

ho free-stream static pressure (approach section, fig. 2), ft of 

liquid abs 

J mechanical equivalent of heat, 778 (ft) (lb force) /Btu 

Kc.min developed cavitation paramctiT based on minimum cavity 

pressure, (/lo— /*<-,min)/(FoV2g) 

Kv developed cavitation parameter based on free-stream vapor 
pressure, (ho—h,)/{Vo^/2g) 

k thermal conductivity of saturated liquid, Btu/ (ft) (hr) (°R) 

L latent heat of vaporization, Btu/ (lb mass) 

NPSH net positive suction head, ft of liquid 

T liquid temperature, °R 

AT decrease in local equilibrium temperature due to vaporization, °R 
Vo free-stream velocity (approach section, fig. 2), ft/sec 
VI volume of saturated liquid, cu ft 

Vv volume of saturated vapor, cu ft 

X axial distance from minimum noncavitating pressure location 

(also cavity’s leading edge) (fig. 2) , in. 

Ax length of cavitated region, in. 

a thermal diffusivity. A:/ (p;C(), (sqft)/hr 

Pi saturated liquid density, (lb mass)/(cu ft) 

saturated vapor density, (lb mass) / (cu ft) 

Subscripts 

max maximum 

min minimum 

pred predicted 

ref reference value obtained by experiment 
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Superscripts 

m, n, p exponents (see eq. (4)) 
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DISCUSSION 


J. HORD (National Bureau of Standards) : The author and his co- 
workers are to be commended for extending the 5-factor* theory and 
synthesizing previous efforts (refs. 1, 2, 12, D-1, D-2) to develop a 
practical and effective solution to a complex and difficult problem. The 
“quasi-static” theory upon which the solution is based is amazingly 
simple, and one cannot help but wonder why it works so well (ref. 14) ; 
nevertheless, it must be concluded that the solution, though simplified, 
contains the significant groups of dimensionless parameters. A current 
research program (ref. D-3) is endeavoring to upgrade the thermostatic 
approach outlined in this paper and account for dynamic effects via 
entrainment theory. Test results from this program indicate that Kc, min 
for a zoro-caliber ogive varies with cavity length as suspected by Ruggeri. 

I feel that the author should warn the reader that equation (3) may 
be a rather crude expression for relating A/i„ and 5 when the fluid proper- 
ties vary appreciably with temperature; as evidenced by the nonlinearity 
of the curves in figure 16, this is particularly true for liquid hydrogen. 
Solution of equation (2) in incremental form as outlined by the author in 
reference 8 is more appropriate for general use. Still a more direct approach 
for computing 5-factor is given below: Hollander (ref. 13) showed that 
various 5-factor solutions (refs. 12, D-1, D-2) were equivalent and could 
be represented by a “quasi-static” model; this model features a unit mass 
of liquid confined in an insulated cylinder with a nonleaking piston. When 
the piston is lifted, wo will assume that the fluid undergoes an isentropic 
expansion. The mass-entropy balance may be written 

miSi = m/2S/2-b»n„2S„2 (D-1) 

and conservation of mass requires that 

mi = nif2+m^2 (D-2) 

Combining equations (D-1) and (D-2) and multiplying the result by 
Pn/ Pv 2 , we obtain 
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where 

mi= initial mass of liquid 
OT /2 = mass of saturated liquid at final pressure 
w»i >2 = mass of saturated vapor at final pressure 
Si = specific entropy of saturated liquid at initial pressure 
S /2 = specific entropy of saturated liquid at final pressure 
Sv 2 = specific entropy of saturated vapor at final pressure 
p /2 = density of saturated liquid at final pressure 
p „2 = density of saturated vapor at final pressure 

This result does not depend on summation of incremental steps to ac- 
count for changing fluid properties, vaporization of liquid, recondensation 
of vapor, and so on. The results are as accurate as the fluid P-V-T data; 
charts, tables, etc., may easily be developed from appropriate equations 
of state or tabulated thermodynamic properties. 

The author emphasizes that local thermodynamic equilibrium within 
the cavitated region is a prerequisite for application of the B-factor 
theory. The “quasi-static” model requires only that thermodynamic 
equilibrium prevail during the vaporization process; consequently, only 
the leading edge of the developed cavity must be in thermodynamic 
equilibrium. The central and trailing edges of the cavity, where condensa- 
tion occurs, may not be in stable thermodynamic equilibrium. Recent 
experiments at NBS with the 0.7-scale venturi substantiate metastability 
in the central and aft regions of vaporous hydrogen cavities. The experi- 
mentally determined exponents in equation (5) are based on pressures 
measured within the cavity near the leading edge, where equilibrium 
prevailed in all fluids tested. It appears that the experimental conditions 
under which the exponents were determined are compatible with the 
thermostatic theory, and local equilibrium throughout the cavitated region 
is not required. It also appears plausible that exponents derived in this 
manner could easily account for slight metastabilities (or perhaps un- 
stable equilibria) within vaporous cavities. 

I heartily agree with the author that more data are needed to sub- 
stantiate the exponent for D/Dret in equation (5). For the scale factor 
studied, the diameter-dependence is less than 4 percent — combining the 
last two terms in equation (5) . For geometrically similar cavities, the last 
term in equation (5) goes to unity and the diameter-dependence is still 
less than 8 percent. 

M. L. BILLET and J. W. HOLL (The Pennsylvania State University) : 
The research reported by the author constitutes a valuable contribution 
to the understanding of thermodynamic effects on developed cavitation. 
Mr. Ruggeri is to be complimented for his measurement of the effect 
and its correlation with possible scaling parameters. Also, the corre- 
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lation of data by means of the S-factor theory is very encouraging. Our 
comments are primarily restricted to introducing an alternate approach 
to the B-factor method based on convective heat transfer (ref. D-3). 

From the B-factor theory of Stahl and Stepanoff (ref. 12) , the tem- 
perature depression is given by 


where B is defined as the ratio of vapor volume to liquid volume and must 
be determined in terms of scaling parameters of the system. However, the 
vapor volume is proportional to the product of cavity length and cavity 
thickness, and also the liquid volume is proportional to the cavity length 
and thickness of the liquid layer supplying the heat of vaporization. Thus, 
B must be a function of a maximum cavity diameter divided by the liquid 
thickness for a body of revolution. 

Previous investigators have shown that for water at ordinary tempera- 
tures the ratio of maximum cavity diameter to model diameter is solely a 
function of the cavitation number based on cavity pressure, which in turn 
is a function of the geometry. 

According to the analysis of Gelder, Ruggeri, and Moore (ref. 8), based 
on the analysis of Eisenberg and Pond (ref. D-4), the liquid thickness is 
proportional to \/«^^/To. Assuming that the cavities are geometrically 
similar, the expression for the B-factor for an axisymmetric flow is 

where Vo is the velocity in the uniform stream. It was assumed that C 
was a universal constant and could be evaluated from a reference state. 
Denoting the reference state with the subscript “ref,” it follows that 


C = 


Bref 

{VoD/a)^J{A^y^, 


Employing equation (D-6) in equation (D-5) for C gives 


B = Bref 


, a ) \VoJ \DrJ L(AA/D)ref. 


(D-6) 


(D-7) 


Application of equation (D-7) to venturi results indicated that the actual 
exponents did not agree with theory. Thus, the exponents in equation 
(D-7) were altered empirically, resulting in equation (5) of the paper. 

An alternate approach to correlating data for the thermodynamic 
effect is the entrainment theory, which is similar to that proposed by 
Acosta and Parkin (ref. D-5) and later bj^ Plesset (ref. D-6). In this 
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theory the latent heat requirement is equated to the convective heat 
transfer at the cavity wall. From this theory, the temperature depression 
is given by 


AT = 


Nu A^/D^ \pj \cj 


where the entrainment coefficient, Cq, is defined as 




Qv 


and where 


(D-8) 


(D-9) 


Aw = surface area of the cavit 3 ^ through which the heat transfer occurs 
Q, = volume flow rate of the vapor in the cavity 
Nu = Nusselt number, hAX/K 
h = film coefficient for convective heat transfer 
Pe= Peclet number, V^AX/oc 

In order to apply equation (D-8), it is necessary to determine the three 
dimensionless parameters Cq, Nu, and A^/D^. 

The cavitation number based on cavity pressure. Pc, is defined as 

where Pq is the free-stream pressure. Similarity for developed cavitation 
in the absence of significant viscous or gravity effects is provided if the 
cavitation number is constant. This result has been verified by numerous 
experiments in water at ordinary temperatures. The similarity law was 
applied to a zero-caliber ogive model (ref. D-3), and the entrainment 
coefficient, Cq, was found by creating a similar ventilated cavity and 
measuring the volume flow rate. The area coefficient, A^/JP, was deter- 
mined from photographs taken of the ventilated cavity. The Nusselt 
number was assumed to be of the form 


Nu = aPe'’P/Fr‘' (D-11) 

where Rc, Pr, and P, are the Reynolds, Prandtl, and Froude numbers, 
respectively. The constants, a, b, c, and d in equation (D-11) were deter- 
mined by substituting the measured values of AT, Cq, and AwjPP and 
equation (D-11) into equation (D-8) and then solving for the constants. 

The equation for the temperature depression for the entrainment 
theory applied to the zero-caliber ogives was found to be 


~o "1 r Zy 1 r 1 w 

AT = constant — — - — yji.sjytMb (D-12) 

LpiJ [Cl\ L“J L^lJ L D J 
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The corresponding equation for the B-factor theory applied to the venturis 
as employed in the paper is 


AT = constant 






(D-13) 


The dimensionless constant in equation (D-12) is equal to 0.495 XlO^’, 
whereas the constant in equation (D-13) is dimensional. 

A comparison of equations (D-12) and (D-13) indicates that the 
j5-factor theory applied to the venturis results in ATa{AX/D)'’-Wo°-^, 
whereas the entrainment theory applied to the ogive yields 
ATa(AX/D)® “Fo“'’- Thus the effect of velocity is greater for the venturis, 
but the effect of cavity length to diameter ratio, AX/D, is greater for 
the ogive. 

Referring to equations (D-12) and (D-13), it is seen that ATaD"-^ for 
the R-factor theory applied to the venturis, whereas ATaZ)® *“ for the 
entrainment theory applied to the ogive. Thus the implication is that the 
effect of size is greater for the ogive. However, this has not been verified 
experimentally, since D was not a variable in the ogive test program; 
thus D enters as a constant in equation (D-12) . Neverthidess, in order for 
equation (D-12) to be dimensionally correct, it is necessary that the 
exponent for D be 0.825. It is intended that the size effect will be checked 
in subsequent investigations of the ogives. 

Referring to the fluid property effects in equations (D-12) and (D-13) , 
it is seen that the influence of (pv/pl) and (L/Cl) are the same in both 
cases, but the exponent on a differs. Also, the entrainment theory con- 
tains the additional term (l/;'i)“S which arises from the kinematic 
viscosity appearing in the Reynolds and Prandtl numbers in the assumed 
relation for the Nusselt number. 

In addition to the above-mentioned differences between the results for 
the venturis and ogives is the variation of cavitation number with AX/ D. 
As indicated in the paper, the cavitation number based on minimum 
cavity pressure is independent of AX/D for the venturis, whereas for the 
ogives the cavitation number decreases as AX/D increases (ref. D-3). 

Equation (D-12), which was obtained from tests with a zero-caliber 
ogive in water from 80° to 290° F w'as employed to predict the behavior of 
a zero-caliber ogive in Freon-113 from 75° to 190° F. The comparisons 
between equation (D-12) and the experiment w'ere found to be very good 
(ref. D-3). As indicated in the paper, equation (D-13) has been applied 
to venturi data obtained with several liquids, and the correlations have 
been encouraging. 

The aforementioned comparisons between the results for venturis and 
zero-caliber ogives show that the effects of fluid properties are not in 
complete agreement and that the major differences appear to be in the 
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effects due to velocity, size, and relative cavity length. This implies a 
need to emphasize the investigation of the influence of changes in bound- 
ary shape and size on the thermodynamic effect. 

F. G. HAMMITT (University of Michigan) : I am interested in whether 
or not the NASA program has included any measurements of void fraction 
in the cavity. The appreciable pressure drop measured across the cavity 
length in the venturi tests seems to indicate the presence of considerable 
liquid within this cavity. I believe such measurements could be obtained 
fairly easily by gamma-ray (or X-ray) densitometry, as we once did in 
our laboratory in a cavitating venturi (ref. D-7) . More detailed informa- 
tion on the mechanics of the actual flow could be very useful in improving 
the realism of the model. 

RUGGERI (author) : I wish to thank all the discussers for their timely 
and very worthwhile comments. First, with regard to Mr. Herd’s dis- 
cussion, I completely agree that equation (3) of the paper does not neces- 
sarily provide a good approximation of the vapor-to-liquid-volume ratio 
curves (R-factors), and one should indeed use caution in applying this 
relation in cases where the physical properties of the liquid and its vapor 
change significantly with temperature. The method for determining 
volume ratio, as outlined by Mr. Hord, is more direct and simple and 
provides ultimate accuracy without the need for solving equations in 
incremental form to account for changes in fluid properties as the tem- 
perature drops due to vaporization. On the question of thermodynamic 
equilibrium within the cavitated region, it should be pointed out that 
although the experimentally determined exponents of equation (5) were 
based on conditions in the cavity’s leading-edge region, where equilibrium 
conditions existed, the final correlations presented in figure 10 included 
considerable liquid-hydrogen data which exhibited various degrees of 
metastability in the central and rearward portions of the cavity. No trend 
was noted in the degree of this metastability. Thus, it appears that 
conditions of thermodynamic equilibrium need only exist in the leading- 
edge region to satisfy the present prediction method. This implies that 
vaporization is confined primarily to the cavity’s leading-edge region, with 
the process for the remainder of the cavity being mostly one of con- 
densation. 

Mr. Billet and Dr. Holl are to be commended for their efforts toward 
the development of a different cavitation model via the entrainment 
theory. The NASA studies were based on a “quasi-static” conduction 
model, which admittedly does not properly represent the real flow situa- 
tion. This may be borne out by the fact that the experimentally deter- 
mined exponents varied somewhat from the theoretical values based on a 
pure conduction process. However, it is interesting to note that, with one 
exception, the parameters for both theories are the same. Nevertheless, 
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the entrainment theory, which is based on a dynamic — or convective — 
model, better represents the real flow situation and thus is more sound 
from a theoretical standpoint. From the standpoint of experimentation, 
however, the entrainment theory involves certain problems similar to 
those encountered with the present cavitation model in the determination 
of absolute rather than effective values of parameters such as Cq and h. 
For example, the volume flow rate measured for a ventilated cavity may 
not necessarily equal the volume flow rate of vapor for a similar vaporous 
cavity. 

The studies of Dr. Holl and his coworkers indicate that for zero-caliber 
ogives, the value of Kc,mm changes with eavity length. This is an inter- 
esting finding, since, during the venturi studies, it had been suspected 
that this variance of Kc, min with cavity length might occur for different 
venturis (or body shapes), but it was not verified. Unquestionably, body 
shape (i.e., noncavitating wall pressure distribution) influences the 
thermodynamic effects of cavitation in general, and this effect is probably 
the primary reason for the difference in exponents between the entrain- 
ment and modified iJ-factor theories. The influence of body shape is 
likely to preclude the prediction of absolute values of thermodynamic 
effects for arbitrary bodies and liquids for considerable time to come. 
Much research is needed in this area, but a more immediate need appears 
to be the influence of body scale on the thermodynamic effects of cavita- 
tion. 

With regard to Dr. Hammitt’s comments, no attempt was made to 
measure the void fraction in vaporous cavities. Attempts were made to 
determine whether macroscopic droplets, wall rivulets, etc., were present 
within the cavity, and, except for isolated instances (particularly for 
liquid nitrogen), none was found. The axial pressure gradients within 
vaporous cavities, in my opinion, are due in large part to the velocity 
(dynamic head) of the vapor, which of course implies very high axial 
velocities in the extreme leading-edge region and comparatively low 
velocities in the collapse region. At the time the NASA Lewis program was 
phased out, plans were to make total and static pressure measurements 
within vaporous cavities to prove or disprove this concept. 
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